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Using the Jordan- Wigner transformation and continued fractions we calculate rigorously the 
thermodynamic quantities for the spin-^ transverse Ising chain with periodically varying intersite 
interactions and/or on-site fields. We consider in detail the properties of the chains having a 
period of the transverse field modulation equal to 3. The regularly alternating transverse Ising 
chain exhibits several quantum phase transition points, where the number of transition points for a 
given period of alternation strongly depends on the specific set of the Hamiltonian parameters. The 
critical behavior in most cases is the same as for the uniform chain. However, for certain sets of the 
Hamiltonian parameters the critical behavior may be changed and weak singularities in the ground- 
state quantities appear. Due to the regular alternation of the Hamiltonian parameters the transverse 
Ising chain may exhibit plateau-like steps in the zero-temperature dependence of the transverse 
magnetization vs. transverse field and many-peak temperature profiles of the specific heat. We 
compare the ground-state properties of regularly alternating transverse Ising and transverse XX 
chains and of regularly alternating quantum and classical chains. 

Making use of the corresponding unitary transformations we extend the elaborated approach to 
the study of thermodynamics of regularly alternating spin-| anisotropic XY chains without field. 
We use the exact expression for the ground-state energy of such a chain of period 2 to discuss how 
the exchange interaction anisotropy destroys the spin-Peierls dimerized phase. 
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I. INTRODUCTORY REMARKS 

The spin-i Ising chain in a transverse field (trans- 
verse Ising chain) is known as the simplest model in 
the quantum theory of magnetism. It can be viewed 
as the ID spin-^ anisotropic XY model in a transverse 
(z) field with extremely anisotropic exchange interac- 
tion. By means of the Jordan-Wigner transformation 
it can be reduced to a ID model of nonintcracting spin- 
less fermionsi* 2 * 3 ^. As a result the transverse Ising chain 
appeared to be an easy case^ and a lot of studies on that 
model have emerged up till now. After the properties of 
the basic skeleton model were understood various modi- 
fications were introduced into the model and the effects 
of the introduced changes were examined. For example, 
an analysis of the critical behavior of the chain with an 



aperiodic sequence of interactions was performed in Ref. 
6, an extensive real-space renormalization-group treat- 
ment of the random chain was reported in Ref. |3, a 
renormalization-group study of the aperiodic chain was 
presented in Ref. ]q It should be remarked, however, that 
the simpler case of the regularly inhomogeneous spin-^ 
transverse Ising chain (in which the exchange interactions 
between the nearest sites and/or the on-site transverse 
fields vary regularly along the chain with a finite period 
p) still contains enough not explored properties which 
deserve to be discussed. Moreover, the thermodynamic 
quantities of such a system can be derived rigorously ana- 
lytically exploiting the fermionic representation and con- 
tinued fractions. 

The thermodynamic properties of the regularly alter- 
nating anisotropic XY chain in a transverse field of pe- 
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riod 2 were considered in Refs. liilOl (see also Ref. [ll| 
where a model without field was investigated). The 
elaborated general approach for calculation of thermo- 
dynamic quantities^ becomes rather tedious if p > 2 
and the properties of chains of larger periods of alterna- 
tion were not discussed. Other papers^^*^ are devoted 
to the ID anisotropic XY model on superlattices, which 
can be viewed as particular cases of a regularly alternat- 
ing anisotropic XY chain. In Ref. the transfer matrix 
method was applied to get the excitation spectrum of 
the Hamiltonian, being a quadratic form of creation and 
annihilation Bose or Fermi operators, on a ID superlat- 
tice. (This fermionic system is related to the ID spin-i 
transverse anisotropic XY model on a superlattice.) In 
Ref. El a version of the approach suggested in Ref. 
was applied to superlattices. Considering as an example 
the ground-state dependences of the transverse magne- 
tization vs. transverse field and of the static transverse 
susceptibility vs. transverse field (which were examined 
numerically for an anisotropic XY chain of period 4) the 
authors of Ref. observed that these quantities behave 
differently than for the isotropic XY model. Contrary 
to the case of isotropic XY model, for the anisotropic 
XY model the number of the critical fields at which the 
susceptibility becomes singular strongly depends on the 
specific values of intersite interaction parameters. The 
quantum critical points in the anisotropic XY chains in 
a transverse field with periodically varying intersite inter- 
actions (having periods 2 and 3) were determined using 
the transfer matrix method in Ref. It was found 

that for periodic chains the number of quantum phase 
transition points may increase and its actual value de- 
pends not only on the period of modulation but also on 
the strengths of anisotropy and modulation of exchange 
interactions. Let us also mention here a paper discussing 
the energy gap vanishing in the dimerized (i.e., period 
2) anisotropic XY chain without fieldi^ and a recent 
paper— which contains such an analysis for a nonzero 
transverse field. 

In the present paper we have obtained a number of 
new results for regularly alternating spin-^ anisotropic 
XY chains exploiting a systematic method for the calcu- 
lation of the thermodynamic quantities not used in the 
references cited above. This approach is based on ex- 
ploiting continued fractions^ and seems to be a natu- 
ral and convenient language for describing regularly al- 
ternating chains (Section [H]). Considering the chains of 
period 3 we examine the generic effects induced by reg- 
ular alternation. We discuss, in particular, the effect of 



regular alternation on the energy gap (Sections [H] and 
llllf) . the zero-temperature dependences of the transverse 
magnetization vs. transverse field and of the static trans- 
verse susceptibility vs. transverse field, and the tem- 
perature dependence of the specific heat (Section lllljl. 
These rigorous analytical results completed by numeri- 
cal calculations of the spin correlation functions (Section 
IIIip demonstrate the effect of the regular alternation of 
Hamiltonian parameters on the quantum phase transi- 
tion inherent in the spin-i transverse Ising chain. Al- 
though in most cases the critical behavior remains like 
in the uniform chain case, for certain sets of the Hamil- 
tonian parameters weak singularities in the ground-state 
quantities may appear. We compare the results for the 
transverse Ising chains with the corresponding ones for 
the isotropic XY chains in a transverse field (transverse 
XX chains); moreover, we also compare the ground-state 
properties of the quantum and classical regularly alter- 
nating transverse Ising/XX chains fSection llII|) . 

The results obtained exploiting the continued fraction 
approach can be used to examine the thermodynamics 
of the regularly alternating spin-^ anisotropic XY chain 
without field since the latter model is related to a system 
of two spin-i transverse Ising chains through certain uni- 
tary transformations. We use the exact expression for the 
ground-state energy of the anisotropic XY chain without 
field of period 2 to demonstrate the effects of anisotropy 
of exchange interaction on the spin-Peierls dimerization 
(Section US- 

We end up this Section introducing notations and mak- 
ing some symmetry remarks. We consider N — * oo spins 
i on a ring governed by the Hamiltonian 

N N 
n—1 n— 1 

(with s N+1 = sf (In = 0) for periodic (open) bound- 
ary conditions). Here I n is the (Ising) exchange in- 
teraction between the nearest sites n and n + 1 and 
f2„ is the transverse field at the site n. We assume 
that these quantities vary regularly along the chain with 
period p, i.e., the sequence of parameters in is 
CI1I1Q2I2 ■ ■ ■ f2 p /pf^ 1/1^2-^2 ■ ■ ■ Qplp ■ ■ ■ ■ Our goal is to 
examine the thermodynamic properties of the spin model 

Let us extend the "duality" transformationiiS to the 
inhomogeneous case (for the Ising chain in a random 
transverse field such a transformation was discussed in 
Refs. It can be easily proved that the partition 

function Z — Trexp (—(3H) for two sequences of param- 
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and . . . /„_if2„ . . . (or . . . I n VL n+ i . . . ) is 
the same. That means that the fields and the interactions 
may be interchanged as f2„ — > I n -i and /„ — > Q„ (or 
&n — * In and /„ — > f2 n +i) remaining the partition func- 
tion unchanged. Really, performing the unitary trans- 
formation [/ = H^L/ ex P {i^SpSp+i) °ne finds that Eq. 

transforms (with the accuracy to the end terms not 
important for the thermodynamics) into 

N N 

UHU+ = J2 + E m ns y n s y n+1 

n—1 ri— 1 

N N 

n—1 n—1 

(to get the second equality we have renumbered the sites 
n — > n — 1 which obviously does not change the ther- 
modynamics). As a result R Z UHU + R Z+ with R z — 
Ylq=i ex P (if s q) ( U P to the end effects) is again the trans- 
verse Ising chain, however, with the exchange interaction 
between the nearest sites n and n+1 being equal to Q n 
(or f2 n +i) and the transverse field at the site n being 
equal to I n -i (or I n ). 

We also recall that the unitary transformation F m = 
2s^ changes the sign of the transverse field at site m 
in the Hamiltonian (|l.lf> . whereas the unitary transfor- 
mation B rn — (2sf) . . . (2sJ„) changes the sign of the ex- 
change interaction between the sites m and m + 1 in 
the Hamiltonian The symmetry remarks permit 

to reduce the range of parameters for the study of the 
thermodynamics of the model. 

Finally, let us extend the relation between the 
anisotropic XY chain without field and the transverse 
Ising chain (see, for example, Refs. E^ll2cll to the inho- 
mogeneous case. Applying the unitary transformation 
V = n^Ti 1 ex P (i 7rs p s p+i) to the Hamiltonian 

N 

F = E (2I» X n+i + 2M**+0 ( L3 ) 

one gets 

N 

VHV+ = E (2^« +2 + I y n s x n+1 ) (1.4) 

n=l 

(with the accuracy to the end terms). This is the Hamil- 
tonian of two independent chains. Performing further 
in Eq. 1)1. 4|) a f rotation of all spins about the y-axis 
one finds that R V VHV + R y+ , R v = JI^li exp (if s v q ) 
(up to the end effects) is the Hamiltonian of two inde- 
pendent transverse Ising chains (in the notations used 



in Eq. where each of ^ sites is defined by the 

sequences of parameters . . . ^+1-^+2-^+3-^+4 ■ ■ ■ and 
. . . ^^+1^+2^+3 ■ ■ ■ • We shall use the discussed rela- 
tion in Section llVl to study the thermodynamic properties 
of regularly alternating anisotropic XY chains without 
field itOl . 

II. CONTINUED FRACTION APPROACH 

To derive the thermodynamic quantities of the spin 
model (11.1(1 we first express the spin Hamiltonian 
in fermionic language by applying the Jordan-Wigner 
transformation^ 2 * 3 -^. As a result we arrive at a model 
of spinless fermions on a ring governed by the Hamilto- 
nian which can be transformed into the diagonal form 

N ( 1\ 
fe=l ^ ' 

after performing a linear canonical transformation. The 
coefficients of the transformation are determined from 
the following equations 1 ' 3 ' 19-21 

ft„_i/ n _i$ fe ,„_i + (nl + i 2 n _ x - AjQ $ fe „ 

+ft„I„$/ 8 , n +i = 0; (2.2) 

+ (n* + In - A fc) *fcn 
+0„+ 1 I n %, n + 1 = 0. (2.3) 

Evidently, we may obtain the thermodynamic quantities 
of the spin model (|1.1|) having the density of states 

1 N 

fe=i 

since due to Eq. (|2.1|l the Hclmholtz free energy per site 
is given by 

/ = -| J™ dEER{E 2 ) In ^2 cosh . (2.5) 

As we shall see below the density of states R(E 2 ) (which 
contains the same information as a set of all A&) is easier 
to calculate than the values of A& or the coefficients $fe„, 

#ftn- 

On the other hand, we may exploit Eqs. I|2.2|l . 12. .'ill 
to obtain the desired density of states R(E 2 ) 1(2.4(1 . We 
note that the three diagonal band set of equations 12. 2|) 
(or 1(2. 3(0 strongly resembles the one describing displace- 
ments of particles in a nonuniform harmonic chain with 
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nearest neighbor interactions and R(E 2 ) 1)2.4)) plays the 
role of the distribution of the squared phonon frequencies 
(for a study of the phonon density of states in a linear 
nonuniform system see, for example, Ref. 122)) . The set of 
equations 1)2.2)1 (or 1)2.3)1 ) can be also viewed as the one 
for determining a wave function of (spinless) electron in 
the ID nonuniform tight-binding model. 

To find the density of states R(E 2 ) from the set of 
equations 12.21) (or 1)2.3)1 ) we use the standard Green func- 
tion approach. Consider, for example, Eq. 12.211 . Let us 
introduce the Green functions G nm = G nm (E 2 ) which 
satisfy the set of equations 



(E 2 — il 2 — I n —i) G nm — ri n _i/ n _iG n _i, 



+ 1,771 



(2.6) 



Knowing the diagonal Green functions G nn = G nn {E 2 ) 
we immediately find the density of states R(E 2 ) 1)2.4)1 
through the relation 

1 N 

R(E 2 ) = V lmG nn (E 2 ± ie), e -» +0. (2.7) 

nN ' 

n=l 

Alternatively, R(E 2 ) can be obtained with the help of 
the Green functions introduced on the basis of the set of 
equations for coefficients ^kn (|2.3|) . The set of equations 
for such Green functions (like Eq. (|2.fi|l ) corresponds 
to the unitary equivalent spin chain (see 1)1.2)1 ) which 
exhibits the same thermodynamic properties. Thus, the 
resulting density of states R(E 2 ) is the same. 

Now we have to calculate the diagonal Green functions 
G nn involved into Eq. I|2.7|l . Let us use the continued 
fraction representation for G nn that follows from 1)2. 6)1 



G r , 



f]2 _ Q2 _ j2 
a si n i n-l 

O 2 T 2 



A 



(2.8) 



E 2 -^-! 2 ^- 



E 2 -n'i 



-t 2 - 



A+ = 



n 2 i 2 



E 2 -n 2 



r- 



E 2 -n 



T 2 

—^72 — 



(Note that the signs of exchange interactions and fields 
are not important for the thermodynamic quantities as 
it was noted above and is explicitly seen from Eq. (|2.8|l .') 
For any finite period of varying f2„ and I n the continued 
fractions in l|2.8|l become periodic (in the limit N — ► oo) 
and can be easily calculated by solving quadratic equa- 
tions. As a result we get rigorous expressions for the 



diagonal Green functions, the density of states 1)2.7)1 and 
the thermodynamic quantities 1)2.5(1 

of the periodically alternating spin chain 1)1. 1)1 . For 
example, on gets for the internal energy e — f + 
for the entropy | = /3 2 §^ or for the specific heat j = 
— /3^jf ■ Assuming that Q n = Q + Afi„ one can also 
obtain the transverse magnetization m z — §^ and the 
static transverse susceptibility \ z = 

Following the procedure described above, for the peri- 
odically alternating chains of period 2 and 3 we find the 
following result for R(E 2 ) 

R(E 2 )={ \M^ 2 ) 2pV ' (2.9) 

0, otherwise, 



where Z p -i(E 2 ) and A2 P (E 2 



nf =1 (E 2 



are 



polynomials of (p — l)th and (2p)th orders, respectively, 
and < a\ < . . . < a 2p are the roots of A2 P (E 2 ). More- 
over, 

Z X {E 2 ) = 2E 2 -Qj-fll- l\ - if, 

p2\ /ir>2r>2 t2 t2 



(E 4 - 



pi + n 2 2 + if + ii) e 



n 2 n 2 + J?7f)t?.io) 



z 2 (e 2 ) = se 4 -2 (nf + nl + n\ + if + 1\ + if) e 2 
+vl\vl 2 2 + n 2 j 2 + i 2 i 2 + njn 2 3 + np 2 + ipl 
+n 2 n 2 + n 2 i 2 + i 2 i 2 , 

~ (e 6 - (n 2 + n 2 2 + n 2 3 + 1 2 + i 2 + z|) e 4 
+ (n\VL 2 2 + ftp 2 + i 2 i 2 + + np 2 + i 2 i 2 
+o 2 o 2 + + ipl) e 2 - njnlnl - /f/f /f X2.11) 

Eqs. 1)2. 9)1 . 1)2. 10)1 . 1)2.11)1 recover the result for the uni- 
form chain if £1„ = f2, /„ = I as it should be. The 
obtained density of states for p = 2 ()2.9I) . 1)2.10)1 can be 
compared with the exact calculation for the anisotropic 
XY chain in a transverse field reported in Ref. EE Such 
a spin chain is represented by noninteracting spinless 
fermions with the energies A±(g) given by Eq. (2.22) 
of that paper. The density of states ()2.4I) has the form 
R(E 2 ) = ±Y^=±J^dq6(E 2 -Al(q)) and for the 
transverse Ising chain of period 2 after a simple inte- 
gration it transforms into 1)2.9)1 . 1)2.10)1 . 

The density of states R(E 2 ) ()2.4)l yields valuable infor- 
mation about the spectral properties of the Hamiltonian 
(11.1)1 . Thus, the gap A in the energy spectrum of the spin 
chain is given by the square root of the smallest root a\ 
of the polynomial A2 P (E 2 ). In Fig. ^ we display the 
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dependence of the energy gap on the transverse field 33 
for some chains of period 2 and 3. The vanishing gap 
indicates quantum phase transition points^. As can be 
seen from the data reported in Fig. ^thc number of such 
quantum phase transition points for a given period of al- 
ternation is strongly parameter-dependent. The chains 
of period 2 may become gapless either at one, two, three, 
or four values of the transverse field, whereas the chains 
of period 3 may become gapless either at one, two, three, 
four, five, or six values of the transverse field depending 
on the specific set of the Hamiltonian parameters. The 
condition for the vanishing gap follows from ^U(O) = 
(|2~TfJ|) (A 6 (0) = HHH) and for the chains of period 2 
(3) it reads il x il 2 = ±l x l 2 (r2i0 2 3 = ±7i7 2 7 3 ). In fact 
we have rederived with the help of continued fractions the 
long known condition for the existence of the zero-energy 
excitations in the inhomogeneous spin-i transverse Ising 
chain^ 3 - which in our notations has the form 

n 1 n 2 ...njv = ±Ji7 2 ...Jiv. (2.12) 

(Notice, that Eq. (6) of Ref. |3 does not contain two 
signs; the minus sign follows from the symmetry argu- 
ments after performing simple rotations of spin axes. It 
is important, as will be seen below, to have two signs in 
l|2.12|l .) Obviously, for periodic chains we have the prod- 
ucts of only p multipliers in the l.h.s. and r.h.s. of Eq. 
(EH1). 

For a chain of period 2 with a uniform transverse field 
Eq. (|2.12JI yields either one critical field il* = if ei- 
ther Ii or 7 2 (or both) equals to zero or two critical fields 
il* = ±y/\lil 2 \. If the transverse field becomes regularly 
varying, ili 2 = il±Ail, Ail > 0, there may be either two 
critical fields il* = ±^/Ail 2 + \hl 2 \ if Ail < y/\hh\, 
or three critical fields il* = (±y/2\hl 2 \, oj if An = 

y/\hh\, or four critical fields il* = ±^ Ail 2 ± |7i7 2 | if 
Ail > ^/|7i7 2 | (see Fig. ^,) . As a result a chain of period 
2 with Ail < y/\lil 2 \, as il varies, exhibits two phases: 
the Ising phase (for |Q| < y/ Ail 2 + \IiI 2 \) and the para- 
magnetic phase (for \il\ > y/ Ail 2 + \IiI 2 \). A chain of 
period 2 with Ail = \/\I\I 2 \, as il varies, also exhibits 
two phases: the Ising phase (for < \il\ < \/2\IiI 2 \) and 
the paramagnetic phase (for \il\ > y'2|7]7 2 |); moreover, 
in the Ising phase at il = il* = the system exhibits a 
weak singularity in the ground-state quantities (see be- 
low). A chain of period 2 with Ail > \f\hT2\, as il 
varies, exhibits three phases: the low-field paramagnetic 
phase (for |0| < \J Ail 2 — |7i7 2 |), the Ising phase (for 
y/Ail 2 - \hh\ < \il\ < y/Ail 2 + |7i7 2 |), and the high- 
field paramagnetic phase (for \il\ > \J Ail 2 + |7i7 2 |). A 



motivation to give such names to different phases follows 
from the behavior of the Ising magnetization m x to be 
discussed below in Section UTTI 

For a chain of period 3 (^2,3 = ^ + Aili 2,3, Afii + 
Ail 2 + AH3 — 0) the critical fields follow from two cubic 
equations 

(O* + Ailx) {il* + Ail 2 ) (il* + Ail 3 ) ± hhh = (2-13) 

each of which may have either one real solution or three 
real solutions. In Fig. [21 we display the regions in AQj- 
Ail 2 plane for the set of parameters of the transverse 
Ising chains with |7i7 2 73| — 1 which yield two (dark re- 
gion), four (gray region), or six (light region) values of the 
critical field. For the set of parameters at the boundary 
between dark and gray (gray and light) regions there are 
three (five) critical fields; for the set of parameters where 
dark, gray and light regions meet there are four critical 
fields. The behavior of the energy gap for all cases can 
be seen in Figs. \flp -QJi. As a result the chain of period 
3 depending on a relation between Af2i, Ail 2 , Ail% may 
exhibit either two phases (the Ising and paramagnetic 
phases), or four phases (two Ising and two paramagnetic 
phases), or six phases (three Ising and three paramag- 
netic phases). Moreover, weak singularities in the Ising 
phases may occur. 

III. THE GROUND-STATE AND 
THERMODYNAMIC PROPERTIES 

A. The ground— state magnetic properties: 
Transverse Ising chain vs. transverse XX chain 

The transverse magnetization and the static transverse 
susceptibility for a regularly alternating transverse Ising 
chain can be obtained using continued fractions as was 
explained in Section ^ Such results for some typical 
chains of period 3 (which roughly correspond to the pa- 
rameters singled out in Figs. \T]p - \$) at zero tempera- 
ture are reported in Fig. [3] Let us compare and contrast 
the results for the magnetic properties of the transverse 
Ising and the transverse XX chains. 

We start from the energy gap. It is known that the 
uniform transverse Ising chain becomes gapless at criti- 
cal field il* — ±|7|. The gap decays linearly while the 
transverse field approaches the critical value, A ~ e, 
e = \il — 0*| — * 0. The transverse XX chain is gap- 
less along the critical line — 17| < il < |/|. The gap opens 
linearly while the value of transverse field exceeds |7|. If 
regular inhomogeneity is introduced into the transverse 
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XX chain the critical line splits into several parts; the 
gaps open linearly as the transverse field runs out the 
critical lines 24 . On the contrary, a regular inhomogene- 
ity introduced into the transverse Ising chain may cither 
only shift the values of critical fields or lead to the ap- 
pearance of new critical points. Moreover, the gap decays 
either linearly, A ~ e, or proportionally to the deviation 
from the critical value squared, A ~ e 2 , as can be seen 
in Fig. n (see also below). 

The energy gap behavior determines the zero- 
temperature transverse magnetization curves for both 
chains. Transverse XX chains exhibit plateaus which 
can be easily understood within the frames of fermionic 
picture. Indeed, a regularly alternating transverse XX 
chain corresponds to a system of free fermions with sev- 
eral energy bands and the transverse field plays the role 
of the chemical potential. Transverse Ising chains do 
not exhibit plateaus, however, being in the paramagnetic 
phases exhibit plateau-like steps (compare the curves in 
Figs. [3t — Ofc and in Figs. QJd -HJl). In the Ising phases 
the transverse magnetization shows a rapid change. In 
the fermionic picture a regularly alternating transverse 
Ising chain again corresponds to a system of free fermions 
with several energy bands, however, the transverse field 
does not play the role of the chemical potential any more. 

The described behavior of the transverse magnetiza- 
tion vs. transverse field is accompanied by the corre- 
sponding peculiarities in the behavior of the static trans- 
verse susceptibility vs. transverse field. Thus, in the 
cases of the transverse XX chain the square-root singu- 
larities indicate the gapless-to-gapped transitions (Figs. 
03 -El). In the case of the transverse Ising chain a lin- 
ear gap decay produces a logarithmic singularity (Figs. 
EH - Gt), whereas for a decay proportional to the squared 
deviation from the critical value the static transverse sus- 
ceptibility does not diverge containing, however, a non- 
analytical contribution which causes a logarithmic sin- 
gularity of its second derivative (short-dashed and long- 
dashed-dotted curves in Figs. Ofi -|2F). 

To end up, we emphasize that for the regularly al- 
ternating transverse XX chains the number of pecu- 
liarities (e.g., in the dependence \ z vs - ^) depends 
only on the period of alternation and equals 2p. This 
is not the case for the regularly alternating transverse 
Ising chains: the number of peculiarities cannot exceed 
2p but may be smaller; the actual number of peculiar- 
ities and their type essentially depends on the specific 
set of the Hamiltonian parameters. Let us also under- 
line a similarity of these results with what has been 



found for the anisotropic/isotropic XY models on ID 
superlattice o 1 ^ 14 ' 2 ^ . 

B. The ground— state magnetic properties: 
Quantum chain vs. classical chain 

To demonstrate the role of quantum effects in the 
zero-temperature magnetization processes we consider 
the classical counterparts of regularly alternating trans- 
verse Ising and transverse XX chains (some calculations 
of the thermodynamic quantities of the uniform classical 
spin chains can be found in Ref. 26). The classical spin 
model consists of classical spins (vectors) s = (s, 0, <j>) 
(0 < < 7r and < <fi < 2ir are the spherical coordinates 
of the spin) on a ring which interact with each other and 
an external field and are governed either by the Hamil- 
tonian 

N 

H = fl n S COS n 

n=l 

N 

+ 2^ ZIjiS 2 sin 9 n sin n +i cos <f> n cos 4> n +i (3.1) 

n=l 

(transverse Ising chain) or by the Hamiltonian 

N 

H = fl n s cos n 

n=l 

N 

+ 2/„s 2 sin0„ sin0 n+ i cos {<j> n - 4>n+x) (3-2) 

n=l 

(transverse XX chain). In Eqs. (|3.1|) . I|3.2|) s is the 

value of the spin which plays only a quantitative role 
(further we put s = i) and the sequence of parame- 
ters for a regularly alternating chain of period p is again 
iif2i/2^2 ■ ■ ■ IpQpIiQ,il2&2 ■ ■ ■ IpQp ■ ■ ■ ■ In what follows 
we restrict ourselves to the case /„ = /, fl n = O + Af2„, 
J2 n Af2„ = which has already been discussed in some 
detail above. Our goal is to examine the effect of regu- 
lar inhomogeneity on the ground-state properties of the 
classical transverse Ising and transverse XX chains. 

Consider at first the transverse Ising chain. One can 
easily construct the ground-state spin configuration and 
the corresponding ground-state energy ansatz. Accord- 
ing to H3.1JI to minimize the ground-state energy one 
should place all spins in xz plane (i.e., <f> n — <j) n +i — 
... = 0(vr) if / < or <\> n = <p n +2 = ... = 0(tt), 
(j> n +i = <f>n+3 = ■■■ = 7r(0) if / > 0). Moreover, the 
angles n are determined to minimize the sum of the con- 
tribution coming from the interaction with the transverse 
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field and of the contribution coming from the intersite in- 
teraction taking into account the period of inhomogene- 
ity. Thus, for the chain of period p, an ansatz for the 
ground-state energy per site reads 



E(e u ...,0 p ) = s ^ 

N p 



^ tin COS 9 n 



2\I\s 
V 



2 P 



sm V n sin V n+ i 



(3.3) 



and the angles 9 n are determined from the set of equa- 
tions 



d E{6 X 



>) 



N 



= 0, n = 1, . 



(3.4) 



Substituting the solution of Eq. (|3.4I) (which yields the 
lowest ground-state energy) into Eq. (|3.3|) we get the 
ground-state energy of the chain. Now the ground-state 
on-site magnetizations are given by m z n — s cos 9 n , = 
s sin n cos (f> n . We can also find the ground-state on-site 
static transverse susceptibility \n — -g^f ■ 

Let us turn to the transverse XX chain 1|3.2I) . In the 
ground-state spin configuration the spin components in 
xy plane are directed arbitrarily but coherently at all 
sites having the values |m^| = ssa\9 n (i.e., <p n — 4> n +i = 
. . . = 4> (</> is an arbitrary angle) if / < or <p n = 4> n +2 = 
... = <i> (0 < <j> < tt), (j) n+ x = 0„ +3 = . . . = (j)+n if I > 0). 
An ansatz for the ground-state energy per site is again 
given by Eq. Ij3.3[l and the angles 9 n are determined from 
Eq. H3.4JI . Moreover, m z n = scos9 n and \n = • 

For the chain of period 1 from (|3.4|) one easily finds 
9 = if uj = jjryy < -1, 9 = arccos (-w) if -1 < uj < 1, 
and 9 = tt if 1 < w. For the chain of period 2 from 
Eq. 13. 4|) in addition to four obvious solutions cos 2 9i = 
cos 2 92 — 1 one gets one more solution 



cos 2 9i = (lu + 5) 



cos 2 2 = (w - sy 



2 l + (u-S)" 



[OJ - 



sy 

1 + (oj + 5) 2 



5 = 



An 
47uT 



(3.5) 



if |w 2 — 8 2 \ < 1. For the chain of period 3 Eq. (|3.4|l has 
again obvious solutions cos 2 9\ = cos 2 9 2 = cos 2 9 3 = 1 ; 
another solution existing at a certain range of the trans- 
verse field can be found numerically (see Figs. -01). 
The described analytical calculations reproduce the re- 
sults obtained earlier numerically for some chains of pe- 
riods 2 and 3 (dashed curves in Figs. 8a, 8b of Ref. 



In Figs. QJi -0J we display the obtained dependences 
of the ground-state magnetizations m z , m x and static 
transverse susceptibility \ z on the transverse field for 
several classical transverse Ismg/XX chains of period 3 
(the results for corresponding quantum chains are shown 
in Fig. OH ■ In contrast to the quantum case, the ground- 
state static transverse susceptibility \ z °f the classi- 
cal chains remains always finite as the transverse field 

varies and hence the classical chains do not exhibit 
any ground-state phase transitions driven by Q. How- 
ever, a regularly alternating classical chain similarly to 
its quantum (XX) counterpart may exhibit plateaus in 
the ground-state dependence transverse magnetization 
m z vs. transverse field fl (compare long-dashed-dotted 
curves in Fig. QJd (classical chain) and Fig. |3h (quan- 
tum chain) which have a plateau — m z = i). Obviously, 
as m z remains constant with varying £1, the static trans- 
verse susceptibility is zero (long-dashed-dotted curves in 
Fig. ^3). Moreover, (m^), n = 1,2, 3 in this region 
is zero (Fig. 0Jr) and the stable ground-state spin con- 
figuration is 9 n = 9 n +i = 7r, 9 n+ 2 = (see Fig. 0]). The 
Ising magnetization m x decays as the system runs out the 
Ising phase according to the power-law, m x ~ |f2* — Q\@, 
with = |. 



C. Quantum phase transitions 

Let us have a closer look at the quantum phase transi- 
tions in regularly alternating transverse Ising chains dis- 
cussing in some detail the critical behavior. For such a 
chain of period p the quantum phase transition points 
are determined by Eq. (|2.12() . The effects of regular al- 
ternation on the number and position of the quantum 
phase transition points in the cases p — 2 and p = 3 were 
analyzed in Section [H] Although Eq. (|2.12l) was found 
many years agoS^ it was not discussed in the context of 
the quantum phase transition theory. In particular, an 
important question how the gap vanishes as the set of 
parameters becomes critical was not considered in Ref. 
1231 Below we show that two types of critical behavior 
are possible: one as it occurs for the second-order phase 
transition (in Ehrcnfest's sense) and another one as it 
occurs for a weaker singularity (the fourth-order phase 
transition in Ehrenfest's sense). These findings are con- 
firmed by numerical computations of the two-site spin 
correlation functions. 

First we analyze how the gap vanishes as the set of 
parameters becomes critical for the case p = 2 when Q* = 
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±y/AQ 2 ± \hl 2 \- Eq. {2~TU)l yields 



A 2 (Q) 



- n*< 



2 ft 2 + Afr 



(3.6) 



If |fi - Q*\ = e -> and ft* ^ Eq. suggests 
A 2 (ft) - e 2 , i.e., the energy gap vanishes linearly (see 
Fig. nji). A linear decay of the energy gap can be also 
seen in many cases in Figs. \T]p -^1 for chains of period 
3. The linearly vanishing gap corresponds to the square- 
lattice Ising model universality class for critical behavior. 
In particular, owing to such a decay of A the ground- 
state energy per site in the vicinity of O* has the form 



/>oo 

e = - / dEE 2 R(E 2 ) 
Jo 



1 

piT 



AEE 1 



f(E 2 



VE 2 - e 2 



-analytical with respect to e 2 terms. 



(3.7) 



Here the first term is a contribution of the lowest energy 
band and the explicit expression for f(E 2 ) is not impor- 
tant for the analysis of nonanalytical behavior as e — > 0. 
The first term in l|3.7|l is proportional to e 2 In e and as a 
result the zero-temperature dependence of m z and x z on 
il contains the nonanalytical terms (fl — f2*) In |fl — f2*| 
and In |fi — respectively. 

Let us turn to the case p = 2 with Ail = v/|ii J2I when 
we have three critical fields fl* = | ± y / 2|7]72T, 1 and 
consider the behavior of A(f2) in the vicinity of Q* = 0, 



i.e., as ft — > 0. From Eq. 



one finds that A 2 (ft) ~ e 4 



Repeating the calculation of the ground-state energy (see 
Eq. IsIIZJl) for such a decay of A one finds that eo con- 
tains the term e 4 In e and hence the system exhibits the 
fourth-order (in Ehrcnfest's sense) quantum phase tran- 
sition at ft* = which is characterized by a logarith- 
mic divergence of the second derivative of the suscepti- 
bility ■ (For an example of a fourth-order thermal 



ample 

f.E3) 



phase transition see Ref. 27.) For p = 3 the dependence 
A(f2) ~ e 2 (see Figs. -CJi) may occur for the sets of 
parameters at the boundaries between different regions 
in Fig. |2| (e.g., at the points b, c, e). Such systems again 
show the fourth-order quantum phase transition behav- 
ior while approaching the corresponding critical fields. 

To discuss further the quantum phases which oc- 
cur as the transverse field varies we examine the spin 
correlation functions (s"s^ +; ). Unfortunately, we can- 
not obtain the spin correlation functions of a regularly 
alternating transverse Ising chain using the continued 
fraction approach which is restricted to the quantities 



that can be expressed through the density of states 
(12. 4|) . However, the spin correlation functions can be 
determined numerically (see, for example, Ref. Esh for 
rather long chains of several thousand sites. Know- 
ing (s^s" +i ) we can obtain the on-site magnetization 



*co\ n°n+rp/ 



Assuming that (s"s£ 



n+rp/ 

( s n)( s n+r P ) decays as (rp)'' 1 exp ^— if r — > oo we 
can also find the correlation length and the power-law 
exponent 7". In our calculations for p = 3 we consider 
chains with N = 2100, take n = 500, rp = 999 to de- 
termine |m|| and rp = 60, . . . , 360 to determine £ x and 
r y x . Our findings are collected in Fig. [S] To illustrate 
the critical behavior of the order parameter in detail we 
also consider a chain of period 2 with = I/2I = 1, 
fii )2 = fl ± 1 taking N = 2000, 4000, 5400, n = f , 
rp = The zero-temperature dependences |mj| vs. ft 
for this chain are reported in Fig. |SJ 

As can be seen from Figs. 03 |B] the behavior of the 
Ising magnetization m x (which plays the role of the or- 
der parameter) indicates the different phases (Ising phase 
if m x ^ or paramagnetic phase if m x = 0) and the 
phase transitions between them. For a set of param- 
eters which yields weak singularities (i.e., m x = in 
the Ising phase, see Figs. |SJd, OJ, Eti Et>) the finite 
size effects are strong and the finite-chain result for x- 
magnetization tends to zero very slowly with increasing 
N (compare data for different N in Fig. HJJd). For the 
second-order (fourth-order) quantum phase transition 
points the critical behavior is given by m x ~ — ill' 3 
with /3 = I (P = \) (compare Figs. and Eh). The 
appearance/disappearance of the Ising magnetization is 
accompanied by a divergence of the correlation length 
£ x = |n — with v = 1 [y = 2) for the second-order 

(fourth-order) quantum phase transition points (Figs. [S^ 
- 0) • Taking into account the values of the exponent 
characterizing the energy gap behavior we conclude that 
the relaxation time scales like the first power of the corre- 
lation length, i.e., the dynamic exponent z = 1, for both 
the second-order and the fourth-order quantum phase 
transitions. At Q = Q* the xx spin correlation functions 
show power-law decay with the exponent "f x = \ for both 
the second-order and the fourth-order quantum phase 
transitions (Figs. [SJn -Et)- Finally, the results for spin 
correlation functions at special values of the transverse 
field f2, i.e., when one on-site field equals zero, coincide 
with the analytical predictions obtained using the 3-site 
cluster Hamiltonian eigenvectors (for a chain of period 2 
the corresponding calculations are given in Ref. |29J) . For 
example, for the chain of period 3 with I\ = I2 = I3 = 1, 
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1, AQo 



0, AQ 3 

'2 



-1 at ft — — 1 we have 
found \mf\ = 5, |mf| « 0.417, |m|| w 0.331 
(see Fig. whereas for such a chain with Afij = 2, 
AO2 = 0, Afi3 = —2 at f2 = —2 we have found 
\mf \ = I, | m || rj 0.267, |mf| « 0.175 (see Fig. Ef). It 
should be noted that the Ising magnetization at the sites 
with zero transverse fields has its maximal value h . Prob- 
ably the most spectacular feature of the Ising chain with 
regularly alternating transverse field is the reentrant be- 
havior with varying nicely seen in Figs. -HJ. The 
appearance of the paramagnetic phase at intermediate 
values of the transverse field when x-magnetization is 
zero and z-magnetization is almost constant can be as- 
sociated with the following classical picture. Assume, for 
example, p = 2 and Q = 0; then owing to the regularly 
varying on-site transverse fields ±Afi with large Aft all 
on-site magnetizations are directed in indirection in 
the spin space. Naturally, this picture may play only an 
auxiliary role for the considered quantum systems. 

Finally we note that our results are in agreement 
with the scaling relations in the theory of conventional 
(temperature-driven) phase transitions^. Thus, the 
quantum phase transition in dimension d = 1 corre- 
sponds to the thermal phase transition in dimension 
d + Z = 2, the exponent v which characterizes the di- 
vergence of the correlation length £ ~ \T — T c \~ w char- 
acterizes the decay of the energy gap A ~ |f2 — 
the exponent a characterizing the divergence of the spe- 
cific heat c ~ |T — T c \~ a characterizes the divergence of 
the transverse susceptibility x z ~ — f2*|~ Q . Moreover, 
a number of scaling relations (which do not account for 
logarithmic divergences) hold. For example, 2 — a = dv. 
Substituting d = 2, v = 1 one gets a = 0, i.e., only a log- 
arithmic divergence in the dependence \ z vs - whereas 
for more rapidly decaying energy gap when v = 2 one 
finds a — —2, i.e., x z does not diverge at 51* (and only 
its second derivative exhibits a logarithmic peculiarity). 



D. Temperature behavior of the specific heat 

We turn to a discussion of the effects of regular al- 
ternation on the temperature dependence of the specific 
heat. The low-temperature behavior of this quantity is 
determined by the fact whether the system is gapped or 
gapless. Thus, the zero-energy excitations immediately 
produce a linear dependence of the specific heat on tem- 
perature. As a result the low-temperature behavior of 
the specific heat indicates the quantum phase transition 
points that can be seen in Fig. [7] in complete agreement 



with the outcomes which follow from the behavior of the 
xx spin correlation functions shown in Fig. More- 
over, we notice that the regular alternation may produce 
many-peak structure of the temperature profiles of the 
specific heat (Fig. EJ. 



IV. ANISOTROPIC XY CHAIN WITHOUT 
FIELD. SPIN-PEIERLS DIMERIZATION 

As a byproduct of the study of regularly alternat- 
ing transverse Ising chains we obtain the thermody- 
namic quantities of regularly alternating anisotropic 
XY chains without field (|1.3|) . Really, using the 
unitary transformations discussed in the end of Sec- 
tion [I] we can state that the Helmholtz free energy of 
the regularly alternating anisotropic XY chain with- 
out field i|1.3fl defined by a sequence of parameters 
Ifimi%...I%iyifimi%...I%iy... is given by Eq. 
(|2.5|) with R{E 2 ) l|2.7|l and the diagonal Green functions 
involved into Eq. Q2.7[l are determined as follows 



E 2 



A" = 



x n-l 

jV 2 jx 1 
L n-2 M n-1 



T y — A~ — A + 



- (4.1) 



E 2 



n — 3 
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A1 



n-2 



I y2 l 



TV 2 jx 2 
J n-4 1 n-3 



E 2 -I* 



n+1 



e 2 - n 



_ TV * 

1 n+2 



-11+ 



Moreover, we may use the obtained densities of states 
(|2.9|) . (|2.10|) . (|2.11|) to find the thermodynamic quanti- 
ties of some regularly alternating anisotropic XY chains. 
Thus, the anisotropic XY chain of period 2 is unitary 
equivalent to two different transverse Ising chains both 
of period 1 and as a result 



R(E 2 ) 



if A xy (E 2 ) > 



2^ ^A IV (E--y 

0, otherwise 

1 1 -,t A IT?2\ 



A aP {E 2 ) 



E 2 



if A yx {E 2 ) > 0, 



otherwise, 



(4.2) 



E 2 



J l ' 1 2 



(Note, that for the isotropic case If = 1% = Ii, 1% 
i| = hi Eq. H4.2[) yields the result obtained in Ref. |2 
(Eqs. (9) - (11) of that paper); in the anisotropic case 
Eq. H4.2|l agrees with the result reported in Ref. Hfl ) 
The anisotropic XY chain of period 3 after performing 
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the above mentioned unitary transformations is equiv- 
alent to two identical transverse Ising chains of period 
3 and, therefore, R{E 2 ) is given by Eqs. (JUJl, 
after the substitution f^i — * If, I\ — ► If, f^2 — * I3, 
I2 —* If, ^3 — > If ' ^3 ~* ^3 • Let us also note that the 
anisotropic XY chains of period 4 (6) are unitary equiv- 
alent to two different transverse Ising chains of period 2 
(3) and hence after simple substitutions Eqs. (|2.9|1 . (|2.1()(l 
(Eqs. H2.9(l . I|2.11|l ) yield the thermodynamic properties 
of such chains. 

Let us use the ground-state energy per site cq = 
—2 AEE 2 R(E 2 ) of the anisotropic XY chain of pe- 
riod 2 to examine the effects of the exchange interac- 
tion anisotropy on the spin-Peierls dimerization inherent 
in the isotropic XY chain^ii 3 ^. For this purpose we as- 
sume in 63 If = (1 + 8){l + 7), If = (1 + S)(l - 7), 
Jf = (1 - <5)(1 + 7), If = (1 - 5)(1 - 7) where < 8 < 1 
is the dimerization parameter and < 7 < 1 is the ex- 
change interaction anisotropy parameter. We consider 
the total energy per site £(8) and its dependence on 8. 
£(5) consists of the magnetic part eo(S) and the elastic 
part a8 2 . Let us recall that in the isotropic limit 7 = 
the total energy £(S) exhibits a minimum at a nonzero 
value of the dimerization parameter 5* 7^ which is a 
manifestation of lattice instability with respect to spin- 
Peierls dimerization—k In the other limiting case 7=1 
the magnetic energy does not depend on 8 and hence the 
uniform lattice is stable. In Fig. [SJi one can see how the 
behavior of £{S) vs. 8 varies as 7 increases from to 0.4 
for a — 0.5. At 7 = the total energy £(§) exhibits a 
minimum at a nonzero value of dimerization parameter 
8* =/= 0. As 7 increases the dependence remains qualita- 
tively the same with, however, slightly decreasing value of 
8* (see Figs. IHb and|St). At a certain value of anisotropy 
parameter ja an additional minimum at 8 — appears. 
Both minima are separated by a maximum, and the min- 
imum at 8* remains the deeper one. At the value 
of 7s (> 1a) the minima have the same depth and with 
further increase of 7 the minimum at 8 = becomes 
the deeper one. If 7 exceeds 7c (> 7s) the minimum 
for a nonzero dimerization parameter disappears. In Fig. 
|SJd one can see the behavior of 8* as 7 varies from 
to 0.4 for different as (solid curves; the dashed curves 
show the behavior of the maximum in the dependence 
£(5) vs. 5); in Fig one can see the dependence of 
8* on 7 for a = 0.4, 0.5, 0.6. The bold dots in the 
curves in this panel correspond to the characteristic val- 
ues of the anisotropy parameter ja < "/b < 7c discussed 
above. The effect of the anisotropy on the spin-Peierls 



dimerized phase occurs according to the first -order phase 
transition scenario. The corresponding phase diagram is 
shown in Fig. [SJi where we indicate the region of stability 
of the dimerized (A) and uniform (C) phases as well as 
the metastable region (regions Bi and B2) where both 
phases may coexist. 

It is interesting to compare the described effects of 
the exchange interaction anisotropy on the spin-Peierls 
dimerized phase with the effects of the transverse field 
on the spin-Peierls dimerized phas o 24 ' 32 . Similarly to the 
anisotropy 7 the transverse field f2 destroys the dimer- 
ized phase according to a first-order phase transition sce- 
nario. However, the value of the dimerization parameter 
8* remains unchanged as f2 increases up to flc- 



V. CONCLUDING REMARKS 

In this work we have analyzed in some detail the 
ground-state and thermodynamic properties of regu- 
larly alternating spin-i transverse Ising chains and 
anisotropic XY chains without field. Due to the Jordan- 
Wigner mapping and the continued fraction approach we 
can calculate the thermodynamic quantities rigorously 
analytically. For certain values of parameters we can also 
calculate the ground-state spin correlation functions. For 
other values of parameters we have calculated the spin 
correlation functions numerically for long chains consist- 
ing of a few thousand sites. We have shown how the 
ground-state properties of regularly alternating classi- 
cal transverse Ising/ AX chains can be examined. The 
main new results obtained are as follows. Firstly, we 
have examined the effects of regular alternation on quan- 
tum phases and quantum phase transitions in the trans- 
verse Ising chain. Owing to regularly alternating param- 
eters the number of quantum phase transition may in- 
crease (but never exceeds 2p where p is the period of 
alternation), the critical behavior remains as in the uni- 
form chain, however, a weaker singularity may also ap- 
pear. Secondly, we have demonstrated how the plateaus 
in the ground-state magnetization curves for the classi- 
cal regularly alternating transverse Ising/AA chains may 
emerge. Thirdly, we have shown how the exchange inter- 
action anisotropy destroys the spin-Peierls dimerization 
inherent in the spin-i isotropic XY chain. The per- 
formed study provides a set of reference results which 
may be useful for understanding more complicated quan- 
tum spin chains. 
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FIG. 1: The energy gap A vs. transverse field Q for transverse 
Ising chains of period 2 (panel a) and 3 (panels b - d) (Q n = 

a + An n , YZ=i An " =o, in = i; a^ 2 = a«i (b), a^ 2 = 

iAfii (c), Afi 2 = -Afii (d)). The bold curves A vs. Q 
correspond to the following values of parameters: Afii — 
0.5, 1, 1.5 (panel a), Afij = 0.5, Afii w 0.630, Afij = 1 
(panel b), Afti = 0.5, Afii « 0.848, Afii = 1.5, Afii £s 
1.921 (panel c), Afii = 1, Afii « 1.375, Afii = 2 (panel d). 
The bold curves in the plane f2-Afh indicate the values of 
parameters which yield the zero-energy gap, A = 0. 
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FIG. 2: Phase diagram of the transverse Ising chain of period 
3 with \hhh\ = 1, ^1,2,3 = 0+Afii,2,3, Afii + Afi 2 +An 3 
J As f2 varies the energy gap vanishes two/four/six times if 
the set of parameters is in the dark/gray/light region. The 
sets of parameters denoted by a, b, c, d, e, f are used below 
to illustrate the dependence on Q of the ground-state Ising 
magnetization (Fig. and the low-temperature specific heat 
(Pig. E| ■ 
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FIG. 3: The ground-state transverse magnetization (a, b, c, 
g, h, i) and static transverse susceptibility (d, e, f, j, k, 1) 
curves for transverse Ising (a - f) and transverse XX (g - 
1) chains of period p = 3. ii,2,3 = 1, fii,2,3 = fi + Afii,2,3, 
Afii + Afi 2 + Afi 3 = 0; a, d, g, j: Afi 2 = Afii, Afii = 
(solid curves), Afii = 0.5 (long-dashed curves), AQi = 0.6 
(short-dashed curves), Afii = 1 (dotted curves); b, e, h, k: 
Afi 2 = -|Afii, Afii = (solid curves), Afii = 0.5 (long- 
dashed curves), Afii = 0.85 (short-dashed curves), Afii = 
1.5 (dotted curves), Afii = 1.9 (long-dashed-dotted curves); 
c, f, i, 1: Afi 2 = -Afii, Afii = (solid curves), Afii = 
1 (long-dashed curves), Afii = 1.35 (short-dashed curves), 
Afii = 2 (dotted curves). 




FIG. 4: The ground-state transverse magnetization (a, b, c), 
Ising magnetization (g, h, i), and static transverse suscepti- 
bility (d, e, f) curves for the classical transverse Ising/XX 
chains of period p = 3. ii,2,3 = — 1, fii,2,3 = fi + Afii,2,3, 
Afii + Afi 2 + Afi 3 = 0; a, d, g: Afi 2 = Afii , Afii = (solid 
curves), Afii — 0.5 (long-dashed curves), Afii = 0.6 (short- 
dashed curves), Afii = 1 (dotted curves); b, e, h: Afi 2 = 
|AOi, Afii = (solid curves), Afii = 0.5 (long-dashed 
curves), Afii = 0.85 (short-dashed curves), Afii = 1.5 (dot- 
ted curves), Afii = 1.9 (long-dashed-dotted curves); c, f, 
i: Afi 2 = -Afii, Afii = (solid curves), Afii = 1 (long- 
dashed curves), Afii = 1.35 (short-dashed curves), Afii = 2 
(dotted curves). We also show the ground-state spin config- 
urations 81, 62, 63 (the corresponding curves are denoted by 
1, 2, 3) for the chains with AQ2 = ^Afii (see panels b, e, 
h) and Afii = 0.85 (j), Afii = 1.5 (k), Afii = 1.9 (1) as fi 
varies. 
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FIG. 5: The ground-state :r-magnetization m x — 
| (\irii + mf | + \m^\) (a, b, c, d, e, f), inverse correlation 
length ^ (g, h, i, j, k, 1), and (m, n, o, p, q, r) for 
the transverse Ising chain of period 3 with |ii| = i2 = 
\h\ = 1, 01,2,3 = n ± Afii, 2 ,3, Afii + An 2 ± AQ 3 = 0, 
Afii = 1, AQ 2 = (fi* w {±1.325}) (a, g, m), Afii = 
1, Afi 2 « 0.292 (fi* « {-1.355,0.678,1.513}) (b, h, n), 
Afii = (f )' « 1.374, Afi 2 = (fi* « {±0.794, ±1.587}) 

(c, i, o), Afii = w 1.374, Afi 2 = -1 (fi* « 

{-1.574,-1.020,-0.348,1.367}) (d, j, p), Afii = 2, Afi 2 « 
-1.575 (O* w {-2.107, -1.874, 0.102, 1.054, 1.772}) (e, k, q), 
Afii = 2, AQ 2 = (0* « {±0.254, ±1.861, ±2.115}) (f, 1, r). 
Connecting curves are guides to the eye. The taken sets of 
parameters are in correspondence with the points a, b, c, d, 
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FIG. 6: The ground-state sublattice :r-magnetizations |mj|, 
j — 1,2 (triangles) and m x = | (\mf | + Imfl) (squares, dia- 
monds, circles) for the transverse Ising chain of period 2 with 
\Ii \ = 1^21 = 1, ^1,2 = fl±l. Connecting curves in panel a are 
guides to the eye. We also report the results for the depen- 
dence m x vs. £1 in the vicinity of the critical points Q* = 
(b) and 0.* — \p2 (c) (squares, diamonds, circles correspond 

to data for N = 2000, 4000, 5400, respectively, solid lines 

1 i 
represent dependences ~ fli (b) and ~ (\/2 — Q) 8 (c)). 




FIG. 7: The low-temperature behavior of the specific heat for 
the transverse Ising chain of period 3 with |7ij = /2I = \h\ = 
1, Qi, 2 ,3 = + Afii,2,3, Afii + AD.2 + AQ 3 = 0, Afii = 1, 
AQ 2 = (a), A«i = 1, A0 2 w 0.292 (b), Afii « 1.374, 
AQ 2 = (c), Afii « 1.374, Afi 2 = -1 (d), Afii = 2, Afi 2 « 
-1.575 (e), Afii = 2, Afi 2 = (f). 
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FIG. 8: The total ground-state energy per site £ (8) vs. 
dimerization parameter S (a = 0.5, from bottom to top 
7 = 0, 0.025, 0.05, 0.1, 0.15, 0.2, 0.4) (a), the dimer- 
ization parameter 8* vs. a (from right to left 7 
J. 0.025, 0.05, 0.1, 0.15, 0.2, 0.4) (b), the dimerization pa- 
rameter 5* vs. 7 (from top to bottom a — 0.4, 0.5, 0.6; the 
meaning of the characteristic values of the anisotropy param- 
eter 7A, 7s, 7c (denoted for a — 0.5) is explained in the main 
text) (c), and the phase diagram in the plane a - 7 (in the 
region A (C) the dimerized (uniform) phase occurs; in the re- 
gions Bi, B2 both phases are possible although in the region 
Bi (B2) the dimerized (uniform) phase is favorable) (d). 



